Defect-induced multicomponent electron scattering in single-walled carbon nanotubes 
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We present a detailed comparison between theoretical predictions on electron scattering processes 
in metallic single- walled carbon nanotubes with defects and experimental data obtained by scanning 
tunneling spectroscopy of Ar + irradiated nanotubes. To this purpose we first develop a formalism 
for studying quantum transport properties of defected nanotubes in presence of source and drain 
contacts and an STM tip. The formalism is based on a field theoretical approach describing low- 
energy electrons. We account for the lack of translational invariance induced by defects within 
the so called extended kp approximation. The theoretical model reproduces the features of the 
particle-in-a-box-like states observed experimentally. Further, the comparison between theoretical 
and experimental Fourier-transformed local density of state maps yields clear signatures for inter- 
and intra- valley electron scattering processes depending on the tube chirality. 

PACS numbers: 73.63Fg,73.22.-f,68.37.Ef 



I. INTRODUCTION 

In the last two decades, carbon nanotubes discovered 
in 1991 1 and in particular Single- Walled Carbon Nan- 
otube (SWNT) that emerged shortly 2,3 after have at- 
tracted a continuously increasing interest from the solid 
states physics community due to their outstanding elec- 
tronic and mechanical properties. 4 9 A single-walled car- 
bon nanotube (SWNT) can be regarded as a strip of 
graphene rolled into a hollow cylinder with diameters 
typically in the order of one nanometer and lengths rang- 
ing from a few microns to a few centimeters. 7 Because of 
the very small diameter, a SWNT is generally considered 
as a quasi-one-dimensional system. SWNTs can be either 
metallic or semiconducting with energy bandgaps varying 
from zero to above 1 eV, depending on their geometrical 
structure. 7 

The low energy spectrum of metallic SWNT comprises 
four branches with a linear dispersion close to the charge 
neutrality point (CNP). Depending on the chirality, the 
corresponding electronic states can be degenerate in the 
orbital momentum. These states are symmetric with 
respect to the energy plane through the CNP imply- 
ing electron- hole symmetry. 10 Other dispersive, energy 
symmetric electronic states are opening at higher energy, 
with a gap > leV. All these electronic states can be ob- 
tained theoretically from a kp method. 6 This approach 
is particularly useful for studying transport properties. 



An ideal metallic SWNT — defectless, charge- neutral 
and infinitely long — can be regarded as a quantum wire 
with four propagating one-dimensional conduction chan- 
nels, i.e. two left- and two right- mo vers. 11 This leads to 
a quantized conductance at zero temperature of 4e 2 //i, 
since each left- or right-moving channel is two-fold de- 
generate in spin. 

In practice, devices with pure and clean SWNTs are 
hardly achievable. Their structure is modified by ran- 
dom potentials arising from adsorbates, defects, inhomo- 
geneous substrate interactions, Schottky barriers at con- 
tacts etc.. These perturbations lead to a finite backscat- 
tering probability and hence an experimentally measured 
conductance below the theoretical maximum of 4e 2 /h. 
The scattering mean-free path £ p of SWNTs has been ex- 
perimentally determined and is of the order of i v ~ 1 /im. 
It is usually longer in metallic SWNTs and shorter in 
semiconducting ones. 7 

SWNTs constitute an ideal platform to study many 
body interactions in one-dimensional conductors. A field 
theoretical approach based on the Tomonaga-Luttinger 
liquid model 12 ' 13 has been proposed a decade ago 14 17 
in order to describe interaction effects in isolated single- 
walled and multi- walled carbon nanotubes. Some exper- 
iments have reported on power-law behavior of the con- 
ductance as function of the applied voltage and/or the 
temperature for SWNTs 18,19 and for metallic junctions 
of crossed SWNTs 20 with exponents in the range of the- 
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oretical predictions. 

Transport experiments on SWNTs cover several dis- 
tinct regimes characterized by the number of external 
electrodes and the quality of the contacts between the 
electrodes and the tube. Source and drain electrodes 
drive an electric current through a contacted SWNT. 
These electrodes can be either weakly or strongly cou- 
pled to the nanotube. In the first case — the so-called 
Coulomb blockade regime — effects of Coulomb charg- 
ing play an important role. 21 Oppositely, for strong cou- 
pling — the so-called Fabry-Perot regime 22 — the trans- 
missivity of the contacts is close to unity 23 25 therefore 
minimizing the effects of Coulomb charging. The use of 
a back gate permits to shift the position of the SWNT 
bands relative to the Fermi level in order to probe the 
conduction channels as a function of energy. In a typical 
strongly coupled two-terminal device, the SWNT differ- 
ential conductance — the derivative of the average cur- 
rent with respect to the applied voltage — shows an os- 
cillatory behavior as a function of the source-drain and 
gate voltages. These oscillations can be easily interpreted 
in terms of a Fabry-Perot electron resonator model: 23 at 
low energy the electrons propagate ballistically between 
the two contacts and the oscillation pattern of the differ- 
ential conductance is the result of coherent interference. 

Detailed investigations of transport and structural 
properties of SWNTs can be achieved by means of scan- 
ning probe techniques as, e.g., scanning tunneling mi- 
croscopy /spectroscopy (STM/STS) or scanning gate mi- 
croscopy (SGM). In both cases, an atomically sharp 
metallic tip is brought in close proximity to the nan- 
otube. In this configuration, by applying a bias voltage 
it is possible to generate a tunneling current between the 
conductive sample and the tip. This current is kept con- 
stant by means of a feedback loop whose output signal 
adjusts the vertical z-position of the tip as a function of 
the x-y position. The recorded z(x,y) signal reflects a 
constant current contour map — it can be interpreted, 
in a first approximation, as the topography of the sam- 
ple. In STS, the tip is stopped at a given position and 
the feedback loop is switched off. The differential con- 
ductance dI/dV{V) is measured by sweeping the bias 
voltage between two values corresponding to two fixed 
energies. The knowledge of the differential conductance 
gives access to the local density of states (LDOS) of the 
sample. 26 

STM/STS techniques allow the study of individual 
SWNTs lying on a conductive substrate 27 30 but also of 
SWNT devices with a source-drain and back-gate geom- 
etry. 31 In SGM, the biased tip — playing the role of a 
local gate — is scanned at fixed height over a SWNT de- 
vice between source and drain contacts. 32 ' 33 The Fermi 
level of the tube is locally shifted due to charge accumu- 
lation near the tip and thus a map of the two-terminal 
transconductance dI/dV g of the nanotube device can be 
established. In all these cases, the tip couples weakly 
to the SWNT. A strong coupling regime can be reached 
with a conductive atomic force microscope (AFM) tip 



brought in contact with the nanotube in order, e.g., 
to perform length-dependent two-terminal conductivity 
measurements . 34 

In this paper, we first investigate theoretically the 
transport properties of strongly coupled SWNT devices 
— perfect source and drain contacts — with a third ter- 
minal, experimentally realized by a conductive tip in the 
weak coupling regime. In absence of translational in- 
variance due to the presence of scattering centers, we de- 
velop a general formalism for studying electron scattering 
within the so-called extended kp model. We account for 
scattering centers such as structural defects, tube ends 
or arbitrary random potentials. The approach leads to a 
general scattering matrix formulation of transport prop- 
erties — current and differential conductance — of im- 
perfect SWNTs. We apply our method to the case of a 
low temperature STM/STS experiment. In this set-up, 
individual finite-length metallic SWNTs dispersed on a 
gold substrate show particle-in-a-box states between con- 
secutive scattering centers. Our approach permits a fun- 
damental understanding of the nature of the observed 
scattering patterns in reciprocal space. 27 

The paper is organized in the following way: in 
Sec. II we give an overview of the electronic properties of 
SWNTs in the kp approximation as well as in the ex- 
tended k p approximation accounting for the presence of 
scattering centers; in Sec. Ill we employ a field theoret- 
ical approach and the extended k p model to develop a 
method suitable to study the quantum transport proper- 
ties of imperfect metallic SWNTs; in Sec. IV we describe 
the experimental setup, the sample preparation, and the 
Fourier transform STS technique used for data analysis; 
in Sec. V we compare the results of theory and experi- 
ment for the case of diverse SWNTs. 



II. THE EXTENDED hp MODEL 

A SWNT can be regarded — from the geometrical 
point of view — as a single-layer graphene strip rolled 
into a cylinder. Its structure is generally indexed by 
its chiral vector defined by the circumferential vec- 
tor that starts and ends on the same lattice site of the 
SWNT (see e.g. Refs. [4,6]). The circumferential vector 
can be expressed as a linear combination of the two basis 
vectors ai = ao(v // 3/2, 1/2) and 3.2 = ao(v // 3/2, —1/2) of 
the honeycomb lattice, where ao = V3a>cc is the lattice 
constant (acc ~ 0.142 nm is the bond length between 
carbon atoms). Thus, = ^iai + ^2 a 2> and the ge- 
ometry of a SWNT is completely defined by the pair of 
integers (711,712) — the chiral indices. The diameter d t 
of a SWNT is given by: 

7 \Ch\ a I o , , 2 /oi\ 

d t 1 — \ n{ + nin 2 + n\ . (2.1) 

The unit cell of a SWNT is defined by the chiral vector 
Ch and the translational vector T perpendicular to C^. 
The translational vector T is the smallest lattice vector 
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that leaves the SWNT lattice invariant when shifted by 
T along the axis and its norm 

T=|T| = ^fV /n ? +nin2+n 2 ( 2 ' 2 ) 

defines the translational period. Here TVr is the great- 
est common divisor of (2ri2 + n\) and (2ni +7^2). The 
nanotube unit cell is thus formed by a cylindrical sec- 
tion of length T and diameter d t . It contains N c = 
4 (n 2 + nm + m 2 ) /TVr carbon atoms. 

To understand the electronic properties of SWNTs, 
a simple way is to start with the band structure of 
graphene. This can be obtained from an effective tight- 
binding model with nearest-neighbor hopping with hop- 
ping integral 70 ~ 2.8 eV. 7 The first Brillouin zone of 
the honeycomb lattice is characterized by the recipro- 
cal lattice vectors bi = (27r/ao)(l/\/3, 1) and b 2 = 
(27r/ao)(l/V / 3, — 1). Only two out of six vertices of 
the first Brillouin zone are non-equivalent — the K 
and K' points. A common choice is usually K = 
(27r/a )(l/v / 3,l/2) and K' = (2ir /a )(l/V3, -1/2). The 
conduction and valence bands formed by 7r-states touch 
at these K and K' points. Due to periodic boundary 
conditions around the circumference of a SWNT, the al- 
lowed wave vectors perpendicular to the tube axis are 
quantized. In contrast, the wave vectors parallel to the 
tube axis remain continuous under the assumption of an 
infinitely long tube. The imposition of periodic boundary 
conditions in transverse direction leads to the following 
constraint on the electronic wavefunction: 

*fe (r + C fc ) = e ifcC " * fc (r) = fc (r) (2.3) 

where r is a vector on the tube surface and k lies in the 
first Brillouin zone of the nanotube. Thus, the electronic 
states are restricted to k vectors that fulfill the condition: 

k c • C h = 27TQ (2.4) 

with q integer. Here kc is the circumferential compo- 
nent of the electron momentum. Plotting the endpoints 
of these allowed vectors for a given SWNT in the re- 
ciprocal space of graphene generates a series of parallel 
and equidistant lines. The distance between the lines 
is Akc = 2/dt. The length, number and orientation of 
these lines depend on the chiral indices (721,712) of the 
SWNT. If we restrict the wave vectors to the first Bril- 
louin zone of the nanotube, its length is found to be 2n/T 
and the number of lines is equal to N c /2. The electronic 
band structure of a specific nanotube is given by the cuts 
of the graphene electronic energy bands along the corre- 
sponding allowed fc-lines: we find a pair of conduction 
(antibonding) and valence (bonding) bands for each fe- 
line (zone folding). Therefore, if the graphene K and K' 
points are crossed by an allowed fc-line, the SWNT is 
metallic, otherwise it is a semiconductor. 

In Fig. 1(a) and (b) we show the case of a (4, 1) SWNT 
in the real and reciprocal space, respectively. Starting 




Figure 1: (Color on-line). Panel (a): the real space lattice 
of a carbon nanotube. The strip rolled into a cylinder is 
displayed as an orange area. The vectors ai and a.2 are a 
basis vector set for the honeycomb lattice, and T are the 
circumferential and the translational vector of the SWNT, 
respectively. Panel (b): The reciprocal space of a carbon 
nanotube is represented by the red line segments. The figure 
correspond to a (4, 1) SWNT — a metallic tube. Here, bi and 
b2 are the reciprocal lattice basis vectors for the honeycomb 
lattice. Panel (c): Energy spectrum for the same SWNT. The 
red-thick lines show the low-energy bands investigated here. 

from the T-point we have a set of 14 parallel lines that 
are touching the vertices of the reciprocal honeycomb lat- 
tice in two points marked with dots — the (4, 1) is then a 
metallic SWNT. 35 In Fig. 1(c) we show the energy spec- 
trum with energies measured relative to the CNP. The 
lines touching the dots representing K and K' of Fig. 1(b) 
are mapped in a set of bands (thick lines) crossing at the 
CNP. Around the CNP these bands have a quasi-linear 
behavior with respect to the momentum. Near the two 
if -points where the energy becomes very small, it is pos- 
sible to perform an effective-mass or kp approximation. 
Here, the wavefunction is factorized in a slowly- varying 
part — containing information about the lattice topology 
— times a plane wave. In this limit, the energy spectrum 
around each of the K points reads 

e±(k) = ±Hv F k (2.5) 

where vp is the effective Fermi velocity — the slope of 
the bands around zero energy. Here k = |&a| is the lon- 
gitudinal component of the momentum relative to the 
if -points. The total longitudinal component of the elec- 
tron momentum is olK ± k (a = ±). The energy spec- 
trum (2.5) presents a four- fold degeneracy: a two- fold 
one due to the presence of two independent if -points, 
and a two-fold one related to the electron spin. This 
degeneracy leads to the fundamental result that the con- 
ductance at zero temperature is quantized and equal to 
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Figure 2: (Color on-line). Reciprocal space for a (4,1) SWNT 
in the case of an impurity breaking the translational invari- 
ance. Here kx and kc are the axial and circumferential axes, 
respectively. The effective size £ s of the scatterer in panel (a) 
is bigger than in panel (b). 



Go = 4e 2 //i for ideal ohmic contacts. 

So far we have considered the properties of a perfect, 
infinitely long SWNT. However, this is the exception, 
in reality we usually are faced with tubes with imper- 
fections, in particular defects in the atomic structure. 
These correspond in some cases to simple vacancies of 
one or more carbon atoms, 7 or lattice distorsions like 
Stone-Wales defects. 36 Adatoms (C, H, N...) also mod- 
ify the chemical structure of SWNTs. 377 In the presence 
of defects, the translational invariance of the SWNTs is 
broken. This implies that the electronic spectrum can- 
not simply be calculated by applying Bloch's theorem as 
done previously. 

In the following we provide an approach tackling the 
problem of quantum transport in SWNTs in the presence 
of defects. First we need to define the important length 
scales of the problem. A defect will distort the lattice 
structure of a SWNT in a region around the defect site 
characterized by a typical length, the effective size £ s of 
the scatterer. This has to be compared with other char- 
acteristic lengths of SWNTs such as the diameter <i t , the 
translational period T, and the atomic lattice spacing a$. 
In addition, defects and impurities give rise to an effective 
scattering potential with a potential strength u s . This 
has to be compared with the effective hopping strength 
7o between carbon atoms. In the following we assume 
that the effective size £ s of the scatterer is smaller than 
the SWNT diameter d t , and that the impurity potential 
u s is smaller than the hopping strength 70. 

Let us start with some considerations concerning the 
smallest of the characteristic length scale of a SWNT: 
the lattice constant ao- In presence of long-range defects 
£ s ^> clq backscattering within the same valley or between 



different valleys is suppressed. This effect — known also 
as Klein tunneling — has been predicted in SWNTs for 
the first time in 1998 by Ando and Nakanishi 39 and has 
recently been observed in SWNT quantum dots. 40 On the 
other hand, if the effective size £ s of the defects is com- 
parable to the lattice constant ao, backscattering within 
the same valley — intra-val\ey scattering — and between 
different valleys — inter- valley scattering — become rel- 
evant. These scattering mechanisms become in general 
exponentially suppressed as the effective size £ s of the 
defect grows. 6 

In order to understand how to go beyond the aforemen- 
tioned standard kp model when dealing with defected 
SWNTs, we have to recall the role played by translational 
invariance in presence of umklapp processes. These pro- 
cesses are characterized by exchanged momenta exceed- 
ing in modulus the length of the reciprocal lattice vectors 
G — the electrons are scattered into a different Brillouin 
zone of the reciprocal lattice. For a translationally invari- 
ant system, the equivalence of wave vectors differing by a 
reciprocal lattice vector allows to reduce these scattering 
events to the first Brillouin zone. 6 

Defects change drastically this physical picture: the 
translational invariance is broken and wave vectors differ- 
ing by reciprocal lattice vectors are no longer equivalent. 
We illustrate this situation in Fig. 2. Here we show the 
honeycomb reciprocal lattice and the parallel lines due 
to the periodic boundary conditions (2.3). In order to 
understand the consequences of the lack of translational 
invariance, we consider two distinct cases: £ s ^> ao in 
panel (a), and ao ~ £ s < T in panel (b). The area of in- 
fluence of defects in the reciprocal space is represented as 
a shadow area of approximative width £~ 1 . The precise 
extent of this area will depend on the details of the shape 
of the Fourier transform of the defect potential. For a sin- 
gle defect, this is in general a decreasing function of the 
distance from the T point. In case (a), scattering events 
induced by the defect produce an exchanged momentum 
smaller than any reciprocal lattice vector. This is usually 
the regime in which backward scattering is suppressed 
— therefore the scattering processes can all be studied 
within the original first Brillouin zone of the translation- 
ally invariant case. In this approximation we can adopt 
the conventional kp picture. However, in case (b) this is 
no longer valid: defects induce scattering processes with 
exchanged momenta exceeding the size of the original 
first Brillouin zone. Hence, we have to account explicitly 
for all processes within a new extended scattering zone 
whose dimensions are given by the inverse effective size 
£~ x of the scatterer. In the extended k p approach we al- 
low for new scattering channels that are associated with 
exchanged momenta larger than the difference between 
the original K points. We refer to this consequences of 
broken translational symmetry as multi-component scat- 
tering. The exchanged momentum 5k of each of these 
scattering processes will be characterized by an axial SkA 
and a circumferential Skc component, respectively (c.f. 
Fig. 4). This approach will be detailed in the following 
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Figure 3: (Color on-line). Outline of the experimental set-up 
and of the Fabry-Perot resonator model. 

section. 



III. MODEL HAMILTONIAN AND 
SCATTERING PROCESSES 

We consider now a metallic SWNT device with ideal 
ohmic contacts to source (S) and drain (D) as well as a 
weak tunneling contact to a third electrode. In our case 
this is an atomically sharp STM tip, henceforth denoted 
as tip (T). The SWNT has a finite length Lq and for the 
x coordinate along it we choose the origin in the middle 
of the nanotube so that the interfaces to the S and D 
electrodes are located at x$ = —Lq/2 and xd = +Xq/2, 
respectively. In the tube we account for the presence of 
several defects characterized by an effective size £ s com- 
parable to the lattice constant ao- The tip is described as 
a semi-infinite non-interacting Fermi gas, and y < de- 
notes the coordinate axis along the tip orthogonal to the 
tube, the origin corresponding to the injection point on 
the tip. The latter is located at position x t with respect 
to the middle of the tube, and electron injection is mod- 
eled by a tunneling Hamiltonian. The total Hamiltonian 
of the system reads 

T~L — ^SWNT + %T + ^tun (3.1) 

where the first term describes the defected SWNT and 
the coupling to the S and D leads. The second term ac- 
counts for the tip, whereas the last one includes SWNT- 
tip tunneling. 

Concerning the SWNT, we focus on the low-energy 
excitations: the energy spectrum is considered as lin- 
ear with respect to the electron momentum around the 
CNP (c.f. Eq. (2.5)). The electron operator ^(x) can 
be decomposed around each of the two Dirac points 
K, K' into right- and left-moving components ^q ?+ (x) 
and ^q ? _(x), respectively. The electron field reads 

tffc) = Yl e i{Q+ak)x ^Qa(x) (3.2) 

Qe{K,K'} 
a=± 

where k denotes the equilibrium Fermi momentum of the 
SWNT with respect to the charge neutrality point Q. 



The Hamiltonian for the SWNT reads 

^SWNT = %kin,SWNT + ^imp + ^/xswnt ■ (3-3) 

In the previous equation the first term 

/CO 
dxa :^Q a (x)d x ^ Qa (x): 

(3.4) 

describes the band linearized spectrum (2.5). The sym- 
bol : : stands for normal ordering with respect to the 
equilibrium ground state. The second term in Eq. (3.3) 
models impurity scattering. In the low-energy approx- 
imation, impurities with an effective size £ s compara- 
ble to the lattice constant ao can be approximated by 
delta-like impurities. In order to fully account for the 
extended kp model introduced in the previous section 
we have to account for multi-component scattering: Each 
of these impurities can either scatter electrons within the 
same valley — intra-valley scattering — or can change 
the valley index of the backscattered electrons — inter- 
valley scattering. In addition, we have to account for 
extended inter-valley scattering processes where the ex- 
changed momenta are larger than the difference between 
the K points. In terms of field operators, the scatter- 
ing Hamiltonian within the extended kp model is of the 
form 

n imp =hv F ^ J2 S Yl A "( A ^> Aa > °) ( 3 - 5 ) 

n K,K' G a,a' 
e -HAK-G + Aak F ] Xn : ^>„)<JwOr„): . 

with AK = K — K' and Aa = a — a' . Here terms with 
a = a' describe forward scattering processes and terms 
with a = —a' describe backscattering. Moreover, for 
K = K' we have intra-valley scattering while terms with 
K 7^ K' describe inter-valley scattering. Finally, G is a 
reciprocal lattice vector and G / describes extended 
inter- valley scattering processes. For a particular case, 
the various scattering processes are sketched in Fig. 4. 
In practice, the coupling constant A n will depend on the 
Fourier transform of the scattering potential of the n- 
th defect. Hence, for a given defect the coefficient A n 
differs for each of the available scattering processes. Also 
the range of G-terms to be summed over depends on the 
spatial extend of the defects. We will allow for several 
defects at positions xi, . . . , xjy- 

The third term in Eq. (3.3) accounts for the source and 
drain electrodes 

/CO 
dx /h S wnt(x)p(x) (3.6) 
-co 

with 

jeVs for x < -x$ ( . 

MSWNt(^) = \ rr r ^ > (3-7) 

I e lor x > xd 

where the term p(x) =: ^(x)^(x) : is the electron density 
fluctuation with respect to the equilibrium value. 
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Figure 4: (Color on-line). Sketch of elastic scattering events 
by a local defect in a chiral (9,0) CNT within the extended kp 
approximation. The parallel (yellow) lines are due to the pe- 
riodic boundary condition along the tube circumference, their 
length is determined by the spatial extend of the defect. The 
(red) circle is the Fermi surface of the valleys at the energy of 
the scattered electrons. Available right [left] moving electron 
states are showed as green [blue] dots. The defect scatters 
electrons between all states indicated as dots. Some typical 
processes are showed explicitly. The green arrow shows an 
intra- valley backscattering event. The blue arrow depicts an 
inter- valley backscattering process. Finally, the black arrow 
indicates an extended inter- valley forward scattering event. 



The Hamiltonian of the tip, the second term in 
Eq. (3.1), reads 



— ^kin,T + • 



(3.E 



Here 



/CO 
dy \c ] (y)d y c(y): (3.9) 
-co 

describes the (linearized) band energy with respect to the 
equilibrium Fermi points ±&t of the tip, and vt denotes 
the Fermi velocity. Notice that the integral runs also over 
the positive ?/-axis, since right and left moving electron 
operators along the physical tip axis y < have been 
unfolded into one chiral (right-moving) operator c(y) de- 
fined on the whole y-axis. The second term in Eq. (3.8) 
describes the bias Vt applied to the tip that affects the 
incoming electrons according to 



/ dy \c ] (y)c(y)\ . 

J — CO 



(3.10) 



Finally, the third term in Eq. (3.1) accounts for the 
SWNT-tip tunneling and reads 

(3.11) 



Htun =hj^/v F V T X 



J2 e ~ KQ+ak)Xt *Qcx( x o) C(0) +h.C. 



where 7 is the tunneling amplitude for electrons, and 
x = xt [y = 0] is the coordinate of the injection point 
along the tube [tip]. 

In the following we will evaluate the differential con- 
ductance in the tip terminal of the described set-up. Ex- 
plicitly, due to the unfolding procedure described above, 
the electron current flowing in the tip at a point y < 
acquires the form 

I T (y,t) = ev T (:c\y,i)c(y,i): - :c\-y,t)c(-y,t):) . 

(3.12) 

Since the system Hamiltonian (3.1) is quadratic in 
the fields ^/Q a (x) and c(y) it admits an exact solu- 
tion. In fact, we can determine the transport prop- 
erties within the Landauer-Biittiker formalism. 22 ' 46 In 
the three-terminal set-up the energy-dependent scatter- 
ing matrix S(E) is a 5 x 5 matrix. It is a function of 
the tip position, the applied tip-SWNT bias, and the 
exchanged momenta 5k in the scattering processes de- 
scribed by the extended k p model. The axial component 
5kA of 5k determines the momentum exchanged by the 
electrons along their direction of motion. On the other 
side, the circumferential component 5kc of 5k affects - 
via the values of A n — the strength of the scattering pro- 
cess. In the following the scattering matrix is expressed 
by the simplified notation 



(3.13) 



where the energy E is measured with respect to the CNP. 
The scattering matrix is obtained by solving the equa- 
tions of motion for the system formed by the SWNT and 
the tip of the STM. The current flowing through the tip 
from the SWNT may be written as 



2e 



dE E(J5) 



(3.14) 



with 



(3.15) 



E(25) = [(|5 51 | 2 + |S 53 | 2 )(/t(£) - fu(E)) 

+(|S 52 | 2 + |S 5 4| 2 )(/t(£)-/s(£))] • 

Here fs(E), /b (.£?), and /t(E) are the Fermi distribution 
functions for the source, drain and tip contacts, respec- 
tively. 

In the low temperature STM/STS experiment de- 
scribed in the next section, the differential conductance is 
deduced from the current flowing through the tip termi- 
nal. In addition, in the experimental set-up the metallic 
SWNT lays on a gold substrate, therefore its potential 
is uniform along the tube length. In our formalism this 
corresponds to the case of chemical potentials of source 
and drain electrodes fixed to the same value, so that 
fs(E) = fr>(E). Given these conditions, the differen- 
tial conductance dl^/dV^ can be readily obtained from 
Eq. (3.14). At zero temperature we have 



dI T 



2e^ 
h 



El* 



'5£\ 



(3.16) 
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This is a function of the tip position and the voltage V? 
applied between the SWNT and the tip. 

We want to add an important remark about the valid- 
ity of the result (3.16) for the differential conductance. 
The strong contact between the metallic SWNT and the 
gold substrate is ruling the phenomena observable exper- 
imentally — this leads to a short lifetime r s of electrons 
injected from the tip into the SWNT. 41,42 Taking into ac- 
count that experiments are performed at 4 K and that the 
applied voltages correspond to even larger energies, ef- 
fects due to interactions with an energy scale smaller than 
hr~ x cannot be resolved in the local density of states. 43 
In addition, the gold substrate is also partially screen- 
ing the electronic Coulomb interaction in the SWNT and 
the tunneling junction resistance is dominating over the 
experimental set-up resistive circuit therefore excluding 
effects as the dynamical Coulomb blockade. 44 Therefore, 
the only effect of interaction that we take explicitly into 
account is a renormalization of the Fermi velocity 45 vp 
in Eq. (3.4). Within this reasonable approximation the 
result (3.16) is valid and will suffice to describe the ex- 
perimental observations. 



IV. EXPERIMENTS 

In the following section, we use our model to study 
the electronic interference patterns observed in the LDOS 
both in real and reciprocal spaces between artificially cre- 
ated scattering centers on metallic SWNTs deposited on 
a conductive substrate. In the next three sections we give 
a brief introduction to the experimental set-up and the 
data analysis technique. 



A. LT-STM/STS setup 

The experiments have been carried out in a commer- 
cial LT-STM set-up (Omicron) operating in ultra high 
vacuum (UHV) with a base pressure below 10 -10 mbar. 
The STM topography images have been recorded in the 
constant current mode. In this set-up we refer to sample 
bias being the potential difference between the grounded 
sample and the tip. We used mechanically cut Pt/Ir tips 
for the topographic scans and the spectroscopy measure- 
ments, where the metallic nature of the tips has been 
regularly checked on the conductive substrate. Differen- 
tial conductance dl/dV spectra — proportional to the 
LDOS in first approximation 26 — were measured using 
the lock-in technique. Under open loop conditions a low 
amplitude a.c. bias (5-15 mV at 600 Hz) is added to the 
sample bias and the resulting a.c. component of the tun- 
neling current is detected. All the measurements have 
been carried out at a temperature of 5 K. All STM topog- 
raphy images presented in this work have been treated 
with the open source WSXM software. 47 



B. Sample preparation 

Gold on glass substrates have been used to support the 
SWNTs for the measurements. In order to obtain clean 
Au (111) monoatomic terraces, the gold surface has been 
prepared by several in situ sputtering and annealing cy- 
cles. The SWNTs provided from Rice University were 
grown by the high pressure CO disproportionation pro- 
cess (HiPco) and purified. 48 To maximize the number of 
individual nanotubes, the raw material was sonicated for 
2-3 h in 1,2-dichloroethane and a droplet of the obtained 
SWNT suspension was then deposited ex situ on the well 
prepared gold surface. Finally, residues from the solvent 
of the suspension were removed in situ by means of a 
short sample annealing at about 390 °C. More details on 
the sample preparation can be found in Ref. [49]. 

Strong electron scattering centers were created by 
two methods: bombardment with medium energy Ar + 
ions inducing vacancy-type defects, 38 ' 52 whereas voltage 
pulses in the STM created cuts of the SWNT at the 
tip position. In the former case, we used a standard 
Leybold™ ion gun designed to be integrated with the 
preparation chamber of the STM setup. The ion gun and 
exposition parameters needed to achieve a defect den- 
sity along the tube axis of about 0.1 nm" 1 have been 
obtained from irradiation tests performed on a freshly 
cleaved graphite (HOPG) substrate, as described in de- 
tails in Ref. [52]. In order to create short SWNT seg- 
ments, we have made use of the technique reported ear- 
lier by Lemay et al. in Ref. [50] where nanotubes were 
cut by applying a local voltage pulse of about \V\ = 5 V 
at constant height (open loop). It is known that the ul- 
trasonication process used to disperse the nanotubes can 
also introduce defects. However, we checked the tubes 
before ion bombardment and found at most one defect 
every 200 nm. Since we aimed at a density of one de- 
fect per 10 nm, not more than one out of twenty defects 
was introduced by the ultrasonication process. The pre- 
cise creation of a scattering defect is in fact not of much 
relevance for the electron scattering processes discussed 
here. 



C. Real and reciprocal space measurements 

Electronic interference patterns induced between con- 
secutive scattering centers are measured by recording a 
series of equidistant STS spectra along the SWNT. In 
a first approximation, these give information about the 
LDOS along the tube axis. 26 Typical dI/dV(x,E(eV)) 
data sets — called di/dV-scans in the following — gener- 
ally consist of 150 dl/dV spectra recorded on topography 
line scans of 300 pts. 

Insights into the electron scattering leading to the in- 
terference patterns are obtained in the reciprocal space 
by means of the Fourier transform STS technique. 
Here, line-by-line zero-padding Fourier transforms of the 
dl / dV -scans are performed in the spatial region of in- 
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Figure 5: (Color on-line), (a) A portion of a metallic SWNT about 40 nm long exposed to a low dose of 1.5 keV Ar + ions 
showing two defects labeled dl and d2. (b) Detail image of the zone delimited by dl and d2. (c) dl /dV-scan recorded along 
the horizontal dashed line in (b). (d) \dI/dV(k, V)\ 2 map resulting from a line-by-line zero padding FFT performed between 
the positions indicated by the two red arrows drawn in (c). (e) Calculated differential conductance for a (14,5) SWNT with 
inter-defect length L — 10.4 nm and asymmetric impurity strengths Ar = 0.35 and Al = 0.15. (f) Corresponding \dl /dV(k, V)\ 2 
map. 



terest — generally between two consecutive scattering 
centers — where the offset amplitude of each line has 
been subtracted. In order to avoid aliasing effects in the 
resulting \dI/dV(k, V)\ 2 data set, we always make sure 
to fulfill the sampling theorem. This requires that the 
highest spectral component of a signal must be smaller 
than 27r/ s /2, where / s is the spatial sampling frequency. 51 
This frequency is defined as / s = l/d s where d s is the 
distance between two consecutive measurement points in 
the dl/dV -scan. For our purpose, we can show that for 
the usual spacing d s = 0.1 — 0.13 nm in a dl /dV -scan, the 
relevant spectral features with k < 25 nm -1 are guaran- 
teed free of aliasing effects up to the maximum analyzed 
length of about 20 nm. 27 

Note that the spatial extent of a dl / dV-scam generally 
acquired within a time frame of about 30 min is always 
observed to be slightly larger (of the order of 5%) than 
the spatial extent of the successive (or previous) topog- 
raphy lines acquired within about 0.5 s, due to a depen- 
dency of the piezoelectric voltage constant on the scan 
velocity. This effect has been systematically corrected 
by a compression of the x-scale of the dl / dV -scan to al- 
low for a consistent comparison between the topographic 



and spatially resolved spectroscopic information on the 
same nanotube. 



V. COMPARISON BETWEEN THEORY AND 
EXPERIMENT 

The first experimental example to which we applied 
our model is illustrated in Fig. 5. Panel (a) shows the 
topography image of a ~ 30 nm long portion of a metallic 
SWNT that has been exposed to a low dose of 1.5 keV 
Ar + ions. 277 A detail image of the zone delimited by 
the pair of defects labeled dl and d2 is displayed in (b). 
Defects induced by medium energy Ar + ions appear typ- 
ically as hillocks with an apparent height ranging from 
0.5 A to 4 A and a lateral extension between 5 A and 
10 A. We have shown in an earlier work 52 , by means 
of a detailed comparison of the spectroscopic signature 
of defect sites with predictions of afr-initio simulations, 
that 1.5 keV Ar + ion bombardment of semiconducting 
SWNTs mainly gives rise to single and double vacancies, 
C-adatoms as well as combinations of these individual 
species. For metallic tubes the same defect structures 
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are expected, but because of prohibitively high compu- 
tational costs, the exact structure of defects in metal- 
lic SWNTs can hardly be determined with high accu- 
racy. However, since the effective size £ s of these defects 
is comparable to the lattice constant ao, they can ap- 
proximately be treated as delta-like impurities with a 
given scattering strength as presented in Sec. III. We 
will show later in this section that this simplified ap- 
proach is accurate enough to interpret our experimental 
data to a large extent. In Panel (c) of Fig. 5 we shows 
the corresponding dl /dV -scan recorded along the hori- 
zontal dashed line in Panel (b). The \dI/dV(k, V)\ 2 map 
resulting from a line-by-line zero padding Fourier trans- 
form performed between the positions indicated by the 
two red arrows drawn in (c) is depicted in (d). The low 
frequency Fourier spots positioned between k = nm" 1 
and k = 4 nm -1 correspond to intra- valley scattering 
processes. 27 They are characterized by an average energy 
separation of AE =160 meV. This is in agreement with 
the value AE =169 meV resulting from the one dimen- 
sional particle-in-a-box model for the linear dispersion of 
a metallic nanotube (E = Hvpk around the inequivalent 
Fermi points). The energy spacing can be obtained by 

* 7r hv F 1.76 Tr ^ 

AE = hv F - = —f- ~ — — eV nm (5.1) 

with L = 10.4 nm and the Fermi velocity v F = 8.5 x 
10 5 m/s. 50 

The more intense dispersive features allowing for in- 
tervalley scattering are visible around k ~ 8 nm -1 and 
at k ~ 16 nm -1 . Further dispersive sets of spots are vis- 
ible around k = 12 nm -1 , k = 21 nm -1 , k = 25 nm -1 
and k = 27 nm -1 . The left-right asymmetry of the in- 
terference pattern in the experimental dl /dV-scan (Fig. 
5c) can be attributed to different scattering strengths of 
the defect sites dl and d2. 27 From atomically resolved 
topography images of the nanotube we can determine a 
chiral angle of « 16° ± 2°. Based on this value of 
the chiral angle and the position of the Fourier spots 
in the \dI/dV(k, V)\ 2 map, we find a good match with 
a (14,5) metallic SWNT (6 = 14.7°). Panel (e) dis- 
plays the calculated differential conductance for a (14,5) 
SWNT with inter-defects length L = 10.4 nm and asym- 
metric impurity strengths Ar = 0.35 and Al = 0.15. 
In general, each of the coefficients \ n (AK, Aa, G) in 
Eq. (3.5) depends on the Fourier transformed scattering 
potential of the defect in question. Here we assume that 
this Fourier transform can adequately be described by a 
simple Gaussian. The interference patterns show char- 
acteristic stripes with increasing curvature when moving 
form the weaker left to the stronger right impurity. 27 All 
the observed Fourier components in (d) are reproduced 
in the calculated \dI/dV(k, V)\ 2 map in (f). The com- 
parison between experiment and theory shows that re- 
garding the positions and slopes of the scattering related 
branches close correspondence is found for the Fourier 
transformed dl/dV- scans. The intensity of the branches 
is, however, not well reproduced by the theory, which 



is the reason for a less pronounced agreement between 
the spatial dl /dV-m&ps. This is mainly due to the fact 
that in reality the scattering strength of a defect will 
depend on energy, momentum and pseudo-spin in a com- 
plex manner, which is not reproduced by the simplified 
delta-like scattering potential of our theoretical model. 
Nevertheless, the agreement between the Fourier trans- 
forms of the theoretical (f) and experimental (d) dl/dV- 
scans clearly proves that the model well reproduces the 
effect of multi-component scattering in metallic SWNTs. 
In the experiment, an asymmetry in the intensity of the 
positive and negative slope branches of the scattering 
process can be observed. In the experiment an asym- 
metry in the intensity of the positive and negative slope 
branches of the scattering process. The dispersion lines 
show a more intense signal for the positive slope branch 
than for the negative one. This asymmetry can be as- 
cribed to a breaking of the particle-hole symmetry due 
to lattice reconstruction that mixes the two carbon sub- 
lattices of the SWNT structure at the defect site. 53 

A second example is illustrated in Fig. 6. Panel (a) dis- 
plays a 23 nm long section of a metallic SWNT irradiated 
with 200 eV Ar + ions showing two defect sites separated 
by L w 17.5 nm. The di/dF-scan shows well-defined 
modes with an energy spacing of about 100 meV as ex- 
pected from the separation of the two defects. In contrast 
to the example above, the interference pattern does not 
show curved stripes. This indicates that in this case the 
scattering strengths of the two defects are similar. There 
is however a discernible symmetry breaking with respect 
to the axis at x = Lo/2, which should be the position 
of a mirror-axis. This points to different scattering phase 
shifts at the two defect sites. Note that our Fabry-Perot 
model does not take the phase shift into account. When 
the tube is cut at both defect sites as shown in panel (b) , 
the mode pattern develops the typical asymmetric stripe 
pattern, which we found to be characteristic for different 
scattering strengths (here with Al > Ar). 27 A further 
cutting pulse close to the left end of the tube results in 
the disappearance of the asymmetry as shown in panel 
(c). It also leads to a weaker resonant modes pattern, 
indicating a reduction of the scattering strength at the 
left end of the tube. The left and right ends have now 
similar scattering strengths. 

This sequence of experiments indicates that even 
though one can expect tube ends to be hard potential 
walls for the electrons, the scattering strengths at tube 
ends can be different. This might result from differ- 
ent effects, e.g., due to different symmetry properties of 
the tube ends or due to a stronger tube-substrate in- 
teraction at the tube end. The latter can be induced 
by the cutting procedure, leading to electron absorp- 
tion instead of scattering. Furthermore, one can notice 
that the \dI/dV(k, V)\ 2 maps corresponding to the three 
cases (a), (b), and (c) show intervalley scattering pat- 
terns with different intensities with some of them almost 
suppressed. To allow a better analysis of these changes, 
we used our model to calculate the differential conduc- 
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Figure 6: (Color on-line). STM topography image and 
corresponding dl /(iV-scan of a metallic SWNT exposed to 
200 eV Ar + ions, (a) Situation with two defect sites. Up- 
per panel shows a raw STM image of an about 30 nm long 
metallic SWNT with three defect sites, lower panel show 
the same data after plane correction and flattening, (b) 
Situation where two cuts have been performed at the de- 
fect positions by means of voltage pulses at « |5| V. (c) A 
third voltage pulse has been performed close to the left end. 
The color scale of the \dI/dV(k, V)\ 2 map is the same for 
the three cases. V s = 0.5 V, I s = 0.3 nA, T = 5.22 K, 
Xres = 0.15 nm. (d) Calculated differential conductance 
with corresponding \dI/dV(k, V)\ 2 map for a (10,4) SWNT 
with inter-defects length L — 20 nm and symmetric impurity 
strengths Ar,l = 0.3. 



tance for a 20 nm long SWNT portion and the corre- 
sponding \dI/dV{k, V)\ 2 map for symmetric scattering 
strengths at both ends (d). The best match was found 
for a metallic (10,4) nanotube. Out of the nine pro- 
jected Fermi cones included in the extended kp model 
[one right/left movers branch with a positive/negative 
slope and a dispersiveless branch corresponding to for- 
ward scattering processes], about six cones are visible in 
the \dI/dV(k, V)\ 2 maps of the three cases presented in 



Fig. 6. In each case, some of the cones are almost sup- 
pressed and some show a higher intensity. Also, in this 
case we observe that some of the presented cones display 
only one dispersive slope. We can interpret these ob- 
servation again in terms of breaking of the particle-hole 
symmetry. In fact in the tube ends several hexagons are 
substituted by pentagons, 4 this results in a direct bond- 
ing of carbon atoms on the same sub-lattice of the SWNT 
structure. 53 



VI. CONCLUSIONS 

We have presented a general method for studying elec- 
tron scattering in metallic SWNTs in presence of defects. 
In particular, we focused on the presence of defects char- 
acterized by an effective size of the order of the SWNT 
lattice constant. In this limit, translational invariance is 
broken implying that the properties of SWNTs cannot 
be obtained from the Bloch theorem. We have shown 
that it is possible to introduce an extended k p model in 
which the momentum exchanged in scattering processes 
can exceed the first Brillouin zone and cover an extended 
region depending on the effective size of defects. Within 
this scheme several scattering channels are available for 
electrons: intra- and inter- valley scattering in addition 
to extended inter-valley scattering. The latter becomes 
relevant because of the lack of translational invariance. 

We have studied the consequences of this kind of de- 
fects on the electronic properties of a metallic SWNT 
contacted by source and drain electrodes. Further, we 
have accounted for a weak coupling to a tip. Based on 
a field theoretical approach describing low-energy elec- 
trons, we have obtained a general relation for the dif- 
ferential conductance between the SWNT and the tip, 
which provides access to the modifications of the local 
density of states induced by defects. We have compared 
the results of our theoretical model for the local density 
of states with the results of experiments where — via the 
STM/STS technique — the same quantity has been mea- 
sured in defected or finite length metallic SWNTs. In the 
experimental set-up defects have been created via Ar + 
ion bombardment, while finite length tubes have been ob- 
tained by means of the pulsed- voltage cutting technique. 
The experimental results clearly show two important fea- 
tures of the theoretical approach. Firstly, the particle-in- 
a-box-like states observed experimentally between con- 
secutive defects are well described by a model with linear 
dispersion for electrons close to the CNP. In accordance 
with the model, the measured level spacing depends only 
on the distance between successive defects. Secondly, the 
extended k p approach fully succeeds in explaining the 
features observed in the Fourier transformed local den- 
sity of states. This includes characteristics coming from 
intra- and inter-valley scattering as well as these associ- 
ated with extended inter-valley scattering. Within this 
model the chiral angle determines uniquely the momenta 
exchanged in the scattering events causing the structure 



11 



of the Fourier transformed local density of states. We 
have shown that the angle obtained from this procedure 
is in agreement with the one measured experimentally by 
atomically resolved topography images of SWNTs. 

In spite of these assets, the model for the impurities 
we have introduced in our theory fails in catching one 
essential feature observed experimentally in the Fourier 
transformed density of states. Specifically, we have rou- 
tinely observed an enhancement of some scattering chan- 
nel with respect to others. These features can be ex- 
plained by a more elaborate atomistic description of de- 
fects accounting for a possible breaking of the particle- 
hole symmetry due to lattice reconstruction at the defect 
sites. 53 
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